Generic dark matter signature for gamma-ray telescopes 
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Abstract 

We describe a characteristic signature of dark matter (DM) annihilation or decay into 
gamma-rays. We show that if the total angular momentum of the initial DM particle(s) van- 
ishes, and helicity suppression operates to prevent annihilation/decay into light fermion pairs, 
then the amplitude for the dominant 3-body final state f + f~"f has a unique form dictated by 
gauge invariance. This amplitude and the corresponding energy spectra hold for annihilation 
of DM Majorana fermions or self-conjugate scalars, and for decay of DM scalars, thus encom- 
passing a variety of possibilities. Within this scenario, we analyze Fermi LAT, PAMELA and 
HESS data, and predict a hint in future Fermi gamma-ray data that portends a striking signal 
at atmospheric Cherenkov telescopes (ACTs). 
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Hopes are high that the long-standing mystery of what comprises the dark matter of our universe 
may be resolved in the coming years. Recent data from the Large Area Telescope (LAT) of the 
Fermi Gamma-ray Space Telescope pQ and from the High Eneregy Stereoscopic System (HESS) [2] 
do not corroborate the excess e + + e~ flux in data from the Advanced Thin Ionization Calorimeter 
(ATIC) [3] and the Polar Patrol Balloon and Balloon borne Electron Telescope with Scintillating 
fibers (PPB-BETS) [I] between 200 and 800 GeV. Also, the LAT 7-ray observations are discrepant 
with Energetic Gamma Ray Experiment Telescope (EGRET) data [5] in that they do not confirm 
the excess at mid-latitudes in the energy range 10 — 50 GeV [6]. Nevertheless, an excess in positrons 
between 10 and 270 GeV in the Payload for Matter Antimatter Exploration and Light-nuclei 
Astrophysics (PAMELA) data [8] may be consistent with the LAT observations. The Fermi 
and PAMELA data may find a common explanation in DM annihilation or decay. If future data 
establish the DM origin, a major breakthrough will be accomplished. 

When a spin-0 DM particle decays, or a pair of identical DM particles of spin-i or spin-0 
annihilate in the static limit, the total angular momentum of the initial state configuration is 
j = 0. Conservation of angular momentum suppresses light fermion pair final states like e + e~ and 
qq without a chirality flip. This suppression disappears if the final state contains an additional 
photon. Following Ref. [9], we refer to the emission of this additional photon in the final state 
as internal bremsstrahlung (IB) to distinguish it from external bremsstrahlung that requires an 
interaction with an external electromagnetic field. IB is comprised of photons radiated from the 
external legs i.e., final state radiation (FSR), and from internal lines i.e., virtual IB (VIbJ^0. In 
this letter, we show that the chirality conserving amplitude for transitions from a j = initial state 
to a e + e~7 final state is given by a unique form governed by QED gauge invariance. It follows that 
for a wide class of DM candidates, a distinct 7— ray signal is expected to accompany the e ± flux. 

We are interested in the chirality preserving amplitude for the final state e(p±) + e{p2) e) ; 
from an initial state of scalar structure: either a decaying DM scalar boson, or a pair of identical 
DM particles with vanishing total angular momentum. The QED gauge invariant amplitude must 
have the form 

M ~ u L (pi)[C(pi,p 2 ) faj^ jk + C(p 2 ,pi) fa/? Pi]vl{P2)^ + (L^R), (1) 

where the form function C becomes a constant if the internal physics is from very short distances. 
Note that the two terms do not interfere in the limit m e — > 0. The operators corresponding to the 



x By a standard abuse of language, FSR is denned to be the leading logarithmic contribution of the photon 
splitting from the external lines and is sometimes erroneously considered gauge-independent. Since VIB is not 
gauge- independent by itself, non-leading logarithmic contributions from external-line-bremsstrahlung need to be 
incorporated to render IB gauge-independent. 

2 The distinction between FSR and VIB is artificial, since in an appropriate gauge, the photons can be thought to 

have been radiated only from the external legs. However, we use this jargon since it is intuitive. 
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Figure 1: Dominant diagrams for the chirality conserving process $ — ► e + e~7, where the initial 
scalar state <3? may be a decaying spin-0 DM particle or a pair of identical DM particles of spin-i 
or spin-0 annihilating with zero total angular momentum. 

amplitude with left chirality are 

where $ is the initial scalar system with mass M$ = 2mr>M for annihilation, and M$ = mr>M for 
decay. The Feynman diagrams corresponding to the process are shown in Fig. [TJ 
The rate for the annihilation process Q^-e + e + jis 

da = (y 2 e) 2 [(l-x 1 ) 2 + (l-x 2 ) 2 ](l-z) 
dx\dz 4-7r 3 m DM (1 — 2x\ — r) 2 (l — 2x2 — r) 2 ' 

where y is the coupling of the intermediate particle of mass the to electrons and the DM parti- 
cle Jf| and r = Am^/M^. The scaling variables x% = 2Ei/M§ for the electron and positron, and 
z = 2E,y/M$ are defined in the static center of mass frame so that x\ + x 2 + z = 2. 

Hard photons arise primarily from VIB from a charged intermediate particle. The photon 
energy distribution is obtained by integrating over x\ G (1 — z, 1) [TO] : 

da __ (y 2 e) 2 1-z / z 2 + (1 + r - z) 2 (1 + r)(l + r - 2z) 1 + r \ 
Wrel dz ~ 327r 3 m 2 DM (1 + r - z) 2 \ Z (1 + r)(l +r - 2i) 1 + r-z " 1 + r - 2z J ' 

If the exchanged particle is much heavier than the DM particle (r — ► oo), the relevant short 
distance physics scale A 3> M$ justifies the use of a dimension-7 operator that is valid for both 
annihilation (with M$ = 2moM) and decay (with M$ = ttidm)'- 

^d v {i> eLl ^ eL )F^ + (L^R) . 

3 As an example, consider a new sector that includes a left-right symmetric electroweak doublet of heavy leptons 
Ll r = (N° , E~)l,r, and a gauge singlet scalar (f>, which is the DM particle. Odd discrete parity is assigned to the 
new particles, and even parity to SM particles. The relevant interaction is £ D U^lLr^), which only applies to the 
left-handed SM lepton doublet l T — (u,e~)L- 
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Z = lEy/MtS, 



z = 2E Y /M<s, 



Figure 2: Photon distributions including VIB compared with the distribution from FSR only 
(dashed curves). The r = 4m|/Mf — > oo case is given by Eq. ([2]). All distributions have unit area 
with the divergent FSR distribution in the left panel cut-off below 0.1 GeV for M<j> = 1 TeV. Note 
the location of the peak in the FSR distribution in the two panels. The distributions in the right 
panel are directly related to Fermi data which are presented as E 2 &j. 

The differential decay distribution is 

- e2Ml -(l-x 3 )[(l-x 1 ) 2 + (l-x 2 ) 2 ] + (L- J R), 



dxidz 512vr 3 A 6 r 



and the total width is 



e 2 M| / 1 1 



+ 



15360tt 3 \A\ A 3 , 

The normalized distributions in z of the prompt photon and in X2 of the positron are 

f| = < 2 > 

~ = 5(3-612 + ^2)12. (3) 

The decay distributions are the same as the normalized annihilation distributions in the short 
distance limit r — » 00. (However, the experimental signatures for decay, and for annihilation with 
r — s- 00, are not identical because the injection flux for decay depends on the DM halo distribution 
as p, while that for annihilation has a p 2 dependence.) The normalized photon distribution is shown 
in Fig. [2j Note that dT/dz increases from low z to peak at z = | and then drops to zero. On the 
other hand, dT /dx2 increases from low X2 of the positron and peaks at the endpoint x^ = 1. 

To reproduce a lifetime of 10 26 s for a DM particle with M$ ~ 2 TeV as suggested by PAMELA, 
Fermi and HESS data, the typical short distance physics scale required is A^ = Ar ~ 10 11 GeV. 



Our study is applicable to the annihilation of Major ana fermions (like neutralinos) and self- 
conjugate scalars, and to decay of DM scalars [11] . Scalar DM annihilation through Higgs exchange 
does not fall within its purview because helicity suppression is not operative; we emphasize that 
the corresponding distributions will be markedly different. Incidentally, in models with spin-1 DM 
like minimal Universal Extra Dimensions (mUED) [12] and Little Higgs with T-parity (LHT) [13J, 
the photon energy distribution (accounting for VIB) is very flat and drops off precipitously at 
m DM [HI- 

We simulate the spectra using GALPROP [15] as described in the appendix of Ref. [16]. The 
default set-up in GALPROP produces p/p spectrum that agrees with PAMELA data [T7] above 
10 GeV for which effects of solar modulation are insignificant. We consider the Einasto [IB] and 
the isothermal [19] DM halo profiles to be representative of mildly cusped and cored profiles, 
respectively. 

In Fig. Owe show that DM annihilation or decay dominantly into e + e~j for M<j> = 1.2 TeV 
easily explains the Fermi LAT, PAMELA and HESS data; larger values of M<j> do not satisfactorily 
reproduce the steeply rising PAMELA positron fraction, and smaller values do not fit the HESS 
data. For DM annihilation, the flux needs to be boosted by about 100 assuming a typical 
thermally averaged annihilation cross section (av) = 3 x 10~ 26 cm 3 s _1 required to produce the 
measured relic abundance. This boost factor is an order of magnitude more than expected from 
N-body simulations [21] . For the decay case the lifetime is about 7 x 10 26 s. Note the significant 
rise in the 7— ray flux at high energies that resembles line emission. The large amplitude of the 
signal distinguishes our class of scenarios from models which produce line emission at loop order. 
Examples with loop-dominated processes include models in which annihilation occurs through s- 
channel Higgs exchange as in several scalar DM models, mUED [22] and LHT [23]. Another 
characteristic is the single large bump; lines from higher order diagrams may appear on the large 
bump with much smaller amplitudes, and may not be resolvable. This is in contrast to models 
that produce multiple lines of roughly equal amplitude [23] , as in theories with two universal extra 
dimensions compactified on a chiral square [25J. Future Fermi data may be able to see a flux 
enhancement with its 10% energy resolution, but will not resolve the bump because it has limited 
sensitivity to photons with energy above 300 GeV. However, ACTs like the CANGAROO III system 
(Collaboration of Australia and Nippon for a Gamma Ray Observatory in the Outback), Major 
Atmospheric Gamma-ray Imaging Cherenkov Telescope (MAGIC), the Very Energetic Radiation 
Imaging Telescope Imaging System (VERITAS), the Cherenkov Telescope Array (CTA), and the 
Advanced Gamma Ray Imaging System (AGIS) which can detect photons in the 50 GeV to 100 TeV 
range, will confirm the bump structure. In Fig. [H we forecast signals for these telescopes in different 
regions of the sky for the two DM halo profiles. We do not consider regions very close to the galactic 
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Figure 3: DM annhilation and decay directly into e + e~7 for M$ = 1.2 TeV. While ACTs will 
clearly see the bump structure in the 7— ray spectrum, Fermi may find a flux enhancement at the 
upper limit of its sensitivity (300 GeV), but will not resolve the bump. Best-fit x 2 values from a 
joint analysis of the PAMELA, Fermi 7— ray, Fermi e^- and HESS datasets which have 7, 18, 26 and 
8 points, respectively are provided in the top-left panel. The number of free parameters is 8, and the 
number of degrees of freedom (dof) is 53, including two energy scale normalizations that account 
for energy calibration uncertainties in the Fermi and HESS datasets. The galactic background 
contribution which was fit for each case is displayed in the same line-type as for the signal. EG is 
the extragalactic 7— ray background that has been estimated up to about 50 GeV [20]. The HESS 
and Fermi error bars have been expanded to approximately include systematic uncertainties (apart 
from the energy scale uncertainties). See the appendix of Ref. [16J for details of the statistical 
analysis. The DM halo follows the Einasto profile. 
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Figure 4: 7— ray signals for DM annihilation with r — > 00 and M$ = 1.2 TeV for the Einasto and 
isothermal profiles in three regions of the sky: inner galactic region, mid-latitudes, and most of the 
sky outside the galactic center region. For both DM profiles, the signals correspond to spectra that 
fit current PAMELA, Fermi and HESS data; in the left panel, very preliminary Fermi data |26j are 
shown for comparison. 

center because the halo profile in the central region is affected by several astrophysical processes, 
that although model-dependent, tend to lower the DM density [27]. A significant bump is seen in 
all regions of the sky for cored and cusped profiles. 

It is known that annihilation of DM Major ana fermions into e is enhanced with a concomitant 
sharp rise in the 7— ray flux close to tudm when electromganetic radiative corrections that relax 
helicity suppression are taken into account [9] . We have shown that the result holds for all scenarios 
in which the total angular momentum of the initial dark matter particle(s) vanishes, and helicity 
suppression operates to prevent annihilation/decay into light fermion pairs, thus applying to scalar 
DM decay, and to static annihilation of identical fermion or scalar DM particles. We also found 
that the chirality preserving amplitude from the initial scalar state to / + /~7 has the unique form 
in Eq. ([1]) dictated by gauge invariance. The experimental signature for a wide class of DM models 
and candidates is unmistakable: a large bump in the high energy 7— ray flux that will be detectable 
by ACTs in any part of the sky. If extant PAMELA, Fermi e and HESS data are to find an 
explanation in DM annihilation or decay, Fermi may see the rising part of the bump. However, if 
M$ is slightly larger than 1.2 TeV, it is entirely possible that Fermi will not see a flux enhancement 
because of its insensitivity to photons above 300 GeV. 
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